3.1370 \(\int (a+b \cos (c+d x))^2 (A+C \cos ^2(c+d x)) \sec ^{\frac{5}{2}}(c+d x) \, dx\)

Optimal. Leaf size=194 \[ \frac{2 \left (a^2 (A+3 C)+b^2 (3 A+C)\right ) \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)} F\left (\left .\frac{1}{2} (c+d x)\right |2\right )}{3 d}-\frac{4 a b (A-C) \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)} E\left (\left .\frac{1}{2} (c+d x)\right |2\right )}{d}+\frac{8 a A b \sin (c+d x) \sqrt{\sec (c+d x)}}{3 d}+\frac{2 A \sin (c+d x) \sec ^{\frac{3}{2}}(c+d x) (a+b \cos (c+d x))^2}{3 d}-\frac{2 b^2 (A-C) \sin (c+d x)}{3 d \sqrt{\sec (c+d x)}} \]

[Out]

(-4*a*b*(A - C)*Sqrt[Cos[c + d*x]]*EllipticE[(c + d*x)/2, 2]*Sqrt[Sec[c + d*x]])/d + (2*(b^2*(3*A + C) + a^2*(
A + 3*C))*Sqrt[Cos[c + d*x]]*EllipticF[(c + d*x)/2, 2]*Sqrt[Sec[c + d*x]])/(3*d) - (2*b^2*(A - C)*Sin[c + d*x]
)/(3*d*Sqrt[Sec[c + d*x]]) + (8*a*A*b*Sqrt[Sec[c + d*x]]*Sin[c + d*x])/(3*d) + (2*A*(a + b*Cos[c + d*x])^2*Sec
[c + d*x]^(3/2)*Sin[c + d*x])/(3*d)

________________________________________________________________________________________

Rubi [A]  time = 0.520299, antiderivative size = 194, normalized size of antiderivative = 1., number of steps used = 7, number of rules used = 7, integrand size = 35, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.2, Rules used = {4221, 3048, 3031, 3023, 2748, 2641, 2639} \[ \frac{2 \left (a^2 (A+3 C)+b^2 (3 A+C)\right ) \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)} F\left (\left .\frac{1}{2} (c+d x)\right |2\right )}{3 d}-\frac{4 a b (A-C) \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)} E\left (\left .\frac{1}{2} (c+d x)\right |2\right )}{d}+\frac{8 a A b \sin (c+d x) \sqrt{\sec (c+d x)}}{3 d}+\frac{2 A \sin (c+d x) \sec ^{\frac{3}{2}}(c+d x) (a+b \cos (c+d x))^2}{3 d}-\frac{2 b^2 (A-C) \sin (c+d x)}{3 d \sqrt{\sec (c+d x)}} \]

Antiderivative was successfully verified.

[In]

Int[(a + b*Cos[c + d*x])^2*(A + C*Cos[c + d*x]^2)*Sec[c + d*x]^(5/2),x]

[Out]

(-4*a*b*(A - C)*Sqrt[Cos[c + d*x]]*EllipticE[(c + d*x)/2, 2]*Sqrt[Sec[c + d*x]])/d + (2*(b^2*(3*A + C) + a^2*(
A + 3*C))*Sqrt[Cos[c + d*x]]*EllipticF[(c + d*x)/2, 2]*Sqrt[Sec[c + d*x]])/(3*d) - (2*b^2*(A - C)*Sin[c + d*x]
)/(3*d*Sqrt[Sec[c + d*x]]) + (8*a*A*b*Sqrt[Sec[c + d*x]]*Sin[c + d*x])/(3*d) + (2*A*(a + b*Cos[c + d*x])^2*Sec
[c + d*x]^(3/2)*Sin[c + d*x])/(3*d)

Rule 4221

Int[(u_)*((c_.)*sec[(a_.) + (b_.)*(x_)])^(m_.), x_Symbol] :> Dist[(c*Sec[a + b*x])^m*(c*Cos[a + b*x])^m, Int[A
ctivateTrig[u]/(c*Cos[a + b*x])^m, x], x] /; FreeQ[{a, b, c, m}, x] &&  !IntegerQ[m] && KnownSineIntegrandQ[u,
x]

Rule 3048

Int[((a_.) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_)*((A_.) + (C_.)*s
in[(e_.) + (f_.)*(x_)]^2), x_Symbol] :> -Simp[((c^2*C + A*d^2)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*(c + d*Sin[
e + f*x])^(n + 1))/(d*f*(n + 1)*(c^2 - d^2)), x] + Dist[1/(d*(n + 1)*(c^2 - d^2)), Int[(a + b*Sin[e + f*x])^(m
- 1)*(c + d*Sin[e + f*x])^(n + 1)*Simp[A*d*(b*d*m + a*c*(n + 1)) + c*C*(b*c*m + a*d*(n + 1)) - (A*d*(a*d*(n +
2) - b*c*(n + 1)) - C*(b*c*d*(n + 1) - a*(c^2 + d^2*(n + 1))))*Sin[e + f*x] - b*(A*d^2*(m + n + 2) + C*(c^2*(
m + 1) + d^2*(n + 1)))*Sin[e + f*x]^2, x], x], x] /; FreeQ[{a, b, c, d, e, f, A, C}, x] && NeQ[b*c - a*d, 0] &
& NeQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0] && GtQ[m, 0] && LtQ[n, -1]

Rule 3031

Int[((a_.) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])*((A_.) + (B_.)*sin[(e
_.) + (f_.)*(x_)] + (C_.)*sin[(e_.) + (f_.)*(x_)]^2), x_Symbol] :> -Simp[((b*c - a*d)*(A*b^2 - a*b*B + a^2*C)*
Cos[e + f*x]*(a + b*Sin[e + f*x])^(m + 1))/(b^2*f*(m + 1)*(a^2 - b^2)), x] - Dist[1/(b^2*(m + 1)*(a^2 - b^2)),
Int[(a + b*Sin[e + f*x])^(m + 1)*Simp[b*(m + 1)*((b*B - a*C)*(b*c - a*d) - A*b*(a*c - b*d)) + (b*B*(a^2*d + b
^2*d*(m + 1) - a*b*c*(m + 2)) + (b*c - a*d)*(A*b^2*(m + 2) + C*(a^2 + b^2*(m + 1))))*Sin[e + f*x] - b*C*d*(m +
1)*(a^2 - b^2)*Sin[e + f*x]^2, x], x], x] /; FreeQ[{a, b, c, d, e, f, A, B, C}, x] && NeQ[b*c - a*d, 0] && Ne
Q[a^2 - b^2, 0] && LtQ[m, -1]

Rule 3023

Int[((a_.) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)] + (C_.)*sin[(e_.) + (
f_.)*(x_)]^2), x_Symbol] :> -Simp[(C*Cos[e + f*x]*(a + b*Sin[e + f*x])^(m + 1))/(b*f*(m + 2)), x] + Dist[1/(b*
(m + 2)), Int[(a + b*Sin[e + f*x])^m*Simp[A*b*(m + 2) + b*C*(m + 1) + (b*B*(m + 2) - a*C)*Sin[e + f*x], x], x]
, x] /; FreeQ[{a, b, e, f, A, B, C, m}, x] &&  !LtQ[m, -1]

Rule 2748

Int[((b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_) + (d_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> Dist[c, Int[(b*S
in[e + f*x])^m, x], x] + Dist[d/b, Int[(b*Sin[e + f*x])^(m + 1), x], x] /; FreeQ[{b, c, d, e, f, m}, x]

Rule 2641

Int[1/Sqrt[sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2*EllipticF[(1*(c - Pi/2 + d*x))/2, 2])/d, x] /; FreeQ
[{c, d}, x]

Rule 2639

Int[Sqrt[sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2*EllipticE[(1*(c - Pi/2 + d*x))/2, 2])/d, x] /; FreeQ[{
c, d}, x]

Rubi steps

\begin{align*} \int (a+b \cos (c+d x))^2 \left (A+C \cos ^2(c+d x)\right ) \sec ^{\frac{5}{2}}(c+d x) \, dx &=\left (\sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)}\right ) \int \frac{(a+b \cos (c+d x))^2 \left (A+C \cos ^2(c+d x)\right )}{\cos ^{\frac{5}{2}}(c+d x)} \, dx\\ &=\frac{2 A (a+b \cos (c+d x))^2 \sec ^{\frac{3}{2}}(c+d x) \sin (c+d x)}{3 d}+\frac{1}{3} \left (2 \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)}\right ) \int \frac{(a+b \cos (c+d x)) \left (2 A b+\frac{1}{2} a (A+3 C) \cos (c+d x)-\frac{3}{2} b (A-C) \cos ^2(c+d x)\right )}{\cos ^{\frac{3}{2}}(c+d x)} \, dx\\ &=\frac{8 a A b \sqrt{\sec (c+d x)} \sin (c+d x)}{3 d}+\frac{2 A (a+b \cos (c+d x))^2 \sec ^{\frac{3}{2}}(c+d x) \sin (c+d x)}{3 d}-\frac{1}{3} \left (4 \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)}\right ) \int \frac{\frac{1}{4} \left (-4 A b^2-a^2 (A+3 C)\right )+\frac{3}{2} a b (A-C) \cos (c+d x)+\frac{3}{4} b^2 (A-C) \cos ^2(c+d x)}{\sqrt{\cos (c+d x)}} \, dx\\ &=-\frac{2 b^2 (A-C) \sin (c+d x)}{3 d \sqrt{\sec (c+d x)}}+\frac{8 a A b \sqrt{\sec (c+d x)} \sin (c+d x)}{3 d}+\frac{2 A (a+b \cos (c+d x))^2 \sec ^{\frac{3}{2}}(c+d x) \sin (c+d x)}{3 d}-\frac{1}{9} \left (8 \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)}\right ) \int \frac{-\frac{3}{8} \left (b^2 (3 A+C)+a^2 (A+3 C)\right )+\frac{9}{4} a b (A-C) \cos (c+d x)}{\sqrt{\cos (c+d x)}} \, dx\\ &=-\frac{2 b^2 (A-C) \sin (c+d x)}{3 d \sqrt{\sec (c+d x)}}+\frac{8 a A b \sqrt{\sec (c+d x)} \sin (c+d x)}{3 d}+\frac{2 A (a+b \cos (c+d x))^2 \sec ^{\frac{3}{2}}(c+d x) \sin (c+d x)}{3 d}-\left (2 a b (A-C) \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)}\right ) \int \sqrt{\cos (c+d x)} \, dx+\frac{1}{3} \left (\left (b^2 (3 A+C)+a^2 (A+3 C)\right ) \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)}\right ) \int \frac{1}{\sqrt{\cos (c+d x)}} \, dx\\ &=-\frac{4 a b (A-C) \sqrt{\cos (c+d x)} E\left (\left .\frac{1}{2} (c+d x)\right |2\right ) \sqrt{\sec (c+d x)}}{d}+\frac{2 \left (b^2 (3 A+C)+a^2 (A+3 C)\right ) \sqrt{\cos (c+d x)} F\left (\left .\frac{1}{2} (c+d x)\right |2\right ) \sqrt{\sec (c+d x)}}{3 d}-\frac{2 b^2 (A-C) \sin (c+d x)}{3 d \sqrt{\sec (c+d x)}}+\frac{8 a A b \sqrt{\sec (c+d x)} \sin (c+d x)}{3 d}+\frac{2 A (a+b \cos (c+d x))^2 \sec ^{\frac{3}{2}}(c+d x) \sin (c+d x)}{3 d}\\ \end{align*}

Mathematica [A]  time = 1.57352, size = 133, normalized size = 0.69 \[ \frac{\sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)} \left (2 \left (a^2 (A+3 C)+b^2 (3 A+C)\right ) F\left (\left .\frac{1}{2} (c+d x)\right |2\right )+\frac{\sin (c+d x) \left (2 a^2 A+12 a A b \cos (c+d x)+b^2 C \cos (2 (c+d x))+b^2 C\right )}{\cos ^{\frac{3}{2}}(c+d x)}-12 a b (A-C) E\left (\left .\frac{1}{2} (c+d x)\right |2\right )\right )}{3 d} \]

Antiderivative was successfully verified.

[In]

Integrate[(a + b*Cos[c + d*x])^2*(A + C*Cos[c + d*x]^2)*Sec[c + d*x]^(5/2),x]

[Out]

(Sqrt[Cos[c + d*x]]*Sqrt[Sec[c + d*x]]*(-12*a*b*(A - C)*EllipticE[(c + d*x)/2, 2] + 2*(b^2*(3*A + C) + a^2*(A
+ 3*C))*EllipticF[(c + d*x)/2, 2] + ((2*a^2*A + b^2*C + 12*a*A*b*Cos[c + d*x] + b^2*C*Cos[2*(c + d*x)])*Sin[c
+ d*x])/Cos[c + d*x]^(3/2)))/(3*d)

________________________________________________________________________________________

Maple [B]  time = 1.561, size = 871, normalized size = 4.5 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((a+b*cos(d*x+c))^2*(A+C*cos(d*x+c)^2)*sec(d*x+c)^(5/2),x)

[Out]

-2/3*(-8*C*(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2)*b^2*cos(1/2*d*x+1/2*c)*sin(1/2*d*x+1/2*c)^6+8*
(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2)*b*(3*A*a+C*b)*sin(1/2*d*x+1/2*c)^4*cos(1/2*d*x+1/2*c)-2*(
-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2)*(A*a^2+6*A*a*b+C*b^2)*sin(1/2*d*x+1/2*c)^2*cos(1/2*d*x+1/2
*c)-2*(2*sin(1/2*d*x+1/2*c)^2-1)^(1/2)*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c
)^2)^(1/2)*(A*EllipticF(cos(1/2*d*x+1/2*c),2^(1/2))*a^2+3*A*EllipticF(cos(1/2*d*x+1/2*c),2^(1/2))*b^2+6*A*Elli
pticE(cos(1/2*d*x+1/2*c),2^(1/2))*a*b+3*C*EllipticF(cos(1/2*d*x+1/2*c),2^(1/2))*a^2+C*EllipticF(cos(1/2*d*x+1/
2*c),2^(1/2))*b^2-6*C*EllipticE(cos(1/2*d*x+1/2*c),2^(1/2))*a*b)*sin(1/2*d*x+1/2*c)^2+A*(-2*sin(1/2*d*x+1/2*c)
^4+sin(1/2*d*x+1/2*c)^2)^(1/2)*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(2*sin(1/2*d*x+1/2*c)^2-1)^(1/2)*EllipticF(cos(1/2
*d*x+1/2*c),2^(1/2))*a^2+3*A*b^2*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(2*sin(1/2*d*x+1/2*c)^2-1)^(1/2)*EllipticF(cos(1
/2*d*x+1/2*c),2^(1/2))*(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2)+6*A*(-2*sin(1/2*d*x+1/2*c)^4+sin(1
/2*d*x+1/2*c)^2)^(1/2)*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(2*sin(1/2*d*x+1/2*c)^2-1)^(1/2)*EllipticE(cos(1/2*d*x+1/2
*c),2^(1/2))*a*b+3*a^2*C*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(2*sin(1/2*d*x+1/2*c)^2-1)^(1/2)*EllipticF(cos(1/2*d*x+1
/2*c),2^(1/2))*(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2)+b^2*C*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(2*sin(
1/2*d*x+1/2*c)^2-1)^(1/2)*EllipticF(cos(1/2*d*x+1/2*c),2^(1/2))*(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)
^(1/2)-6*C*(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2)*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(2*sin(1/2*d*x+1/
2*c)^2-1)^(1/2)*EllipticE(cos(1/2*d*x+1/2*c),2^(1/2))*a*b)/(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2
)/(2*cos(1/2*d*x+1/2*c)^2-1)^(3/2)/sin(1/2*d*x+1/2*c)/d

________________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int{\left (C \cos \left (d x + c\right )^{2} + A\right )}{\left (b \cos \left (d x + c\right ) + a\right )}^{2} \sec \left (d x + c\right )^{\frac{5}{2}}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*cos(d*x+c))^2*(A+C*cos(d*x+c)^2)*sec(d*x+c)^(5/2),x, algorithm="maxima")

[Out]

integrate((C*cos(d*x + c)^2 + A)*(b*cos(d*x + c) + a)^2*sec(d*x + c)^(5/2), x)

________________________________________________________________________________________

Fricas [F]  time = 0., size = 0, normalized size = 0. \begin{align*}{\rm integral}\left ({\left (C b^{2} \cos \left (d x + c\right )^{4} + 2 \, C a b \cos \left (d x + c\right )^{3} + 2 \, A a b \cos \left (d x + c\right ) + A a^{2} +{\left (C a^{2} + A b^{2}\right )} \cos \left (d x + c\right )^{2}\right )} \sec \left (d x + c\right )^{\frac{5}{2}}, x\right ) \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*cos(d*x+c))^2*(A+C*cos(d*x+c)^2)*sec(d*x+c)^(5/2),x, algorithm="fricas")

[Out]

integral((C*b^2*cos(d*x + c)^4 + 2*C*a*b*cos(d*x + c)^3 + 2*A*a*b*cos(d*x + c) + A*a^2 + (C*a^2 + A*b^2)*cos(d
*x + c)^2)*sec(d*x + c)^(5/2), x)

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*cos(d*x+c))**2*(A+C*cos(d*x+c)**2)*sec(d*x+c)**(5/2),x)

[Out]

Timed out

________________________________________________________________________________________

Giac [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int{\left (C \cos \left (d x + c\right )^{2} + A\right )}{\left (b \cos \left (d x + c\right ) + a\right )}^{2} \sec \left (d x + c\right )^{\frac{5}{2}}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+b*cos(d*x+c))^2*(A+C*cos(d*x+c)^2)*sec(d*x+c)^(5/2),x, algorithm="giac")

[Out]

integrate((C*cos(d*x + c)^2 + A)*(b*cos(d*x + c) + a)^2*sec(d*x + c)^(5/2), x)